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Abstract 

Pion Compton scattering is studied in perturbative QCD for real and 
space-like initial photons. Different methods for the convolution of the 
hard amplitude with the pion wave-functions, which have in the past led 
to conflicting results, are compared. 
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Introduction 



The analysis of exclusive processes in QCD has vastly improved in the last few years. 
The role of Sudakov logarithms in elastic scattering has been investigated It has 
been shown || that, like form factors and 77 reactions, large angle photoproduction and 
Compton scattering are not affected in lowest order by Sudakov corrections and that, 
as a consequence, the magnitude and the phase of the scattering amplitudes for these 
processes can be reliably predicted by perturbative QCD (pQCD) at large momentum 
transfer. The applicability of pQCD to hadron form factor has been extended to lower 
energies using an improved version of the factorization formula ||. The sum-rule 
approach has been applied to Compton scattering [|], [5]. 

Nonetheless, the number of reactions for which theoretical predictions exist is still 
rather small. The difficulties are best exemplified by the existence of three published 
calculations for proton Compton scattering, 7p — > jp, in disagreement with each other 

I Si- 

In this paper we consider pion Compton scattering for real and space-like initial 
photons in pQCD. Pion Compton is the simplest reaction in which non trivial-phases 
are present and which is experimentally accessible. Its calculation is conceptually simi- 
lar to the calculation of proton Compton scattering, but much simpler. It is, therefore, 
an ideal laboratory for studying the different methods || |7|, for the convolution of 
the hard amplitude. 

Pion Compton has been considered before for real j5], |9| and time-like photons |§, 
in this paper we extend previous calculations to new kinematical domains and correct 
some errors that are present in |§. 

There are different ways of measuring the 7T7 cross section. In the first instance a 
beam of pions is scattered off of a nucleus. In the forward direction the one-photon- 
exchange mechanism dominates [[UJ and tt'j scattering can be studied. Alternatively 
the interaction can be between an high-energy electron beam and a nucleus, with the 
extraction of quasi-real pions from nuclear matter. 



Calculation 

Computing the amplitude for an exclusive process involves two steps. First of all, an 
analytical expression for the large set of diagrams which is often required needs to be 
derived. This is usually far from trivial and in some cases can only be accomplished 
by some automatic procedure. Second, the result has to be multiplied by the wave 
functions of all external hadrons and then it has to be integrated over the momentum 
fractions of all external quarks. Apart from simple cases like lepton-hadron scattering 
and photon-photon annihilation, for particular values of the quark momenta, one or 
more of the particles propagating inside the various diagrams can go on their mass-shell. 
Mathematically this means that one-dimensional integrals of the form: 

lim f 1 f ^ — dx = P ^ I^-dx-iirfta) (1) 
e-+o Jo x — a + ie Jo x — a 

where < a < 1, have to be computed, and that the result has to be integrated over 
all remaining variables. 



Several methods have been proposed but unfortunately when applied to the reaction 
IP ~ y IP they have led to conflicting results. Therefore we have preferred not to rely on 
only one method of integration. This also allows us to compare the different procedures 
and to study their relative merits. 

We have considered a simpler version of the method used in ||, which we will call 
the subtraction method, the technique suggested in 0, which we will call the KN 
method and the procedure proposed in ||, which will be referred to as the ie method. 

In the first two cases some transformation is applied to the principal part integral 
in the right-hand side of ([!]) in order to avoid the large cancellations involved which 
lead to large numerical instabilities. The imaginary part is then computed by hand, a 
procedure that is rather tedious and that, if not automated, can be a source of errors 
which are almost impossible to detect. 

The subtraction method is based on the following identity: 

pfi^dx = r 1 /(x) - /(a) ^ +/ (q ) p/ 1 ^^ 

Jo x — a Jo x — a Jo x — a 

1 /0) - f{a) . . 1 - a 

dx + f(a) log (2) 

o x — a a 

Both resulting integrals are numerically stable. 

The KN method relies on a change of variable to simplify the evaluation of the real 
part of eq. (0). One can write: 

P f 1 1&L dx = limfJi + J 2 ) (3) 
Jo x — a e-+o 

where 

h = r € dx j 2 = f 1 dx (4) 

JO X — a Ja+e x — a 

Now the following change of variables is performed on J x : 

x T /-W" f{a Vl ) 

Vi = ~ => Ji = / r dyi (5) 

a Jo 2/1 — 1 



and on J 2 : 



where 



Sfc = s , /-, =^ J 2 = / —77^ r~ dy 2 (6) 



d 

The integration variables y%, y 2 satisfy: 



a 2 + (l-2a)x Jo d(l - y) 

9= a(1 ~ ^ d = a + (l-2a)y (7) 



dx , dx 



d yi ^ =1 " dy 2 ^ =1 y ' 

therefore the upper limits in (|^, |6[) approach 1 uniformly. From ([5], |^) one can derive an 
expression for the principal part of the integral in which the limit e — > can be safely 
taken: 



Pf 1 M. dx= Nm+ fi-MW Xdy ( 9 ) 

Jo x — a Jo \y — 1 d(l — y) J 

The last integral is well behaved and can be directly integrated numerically. 

In our experience the subtraction method and the KN method always give the same 
results and can be implemented with comparable effort. 

The ie method is the easiest one to implement automatically on the large set of 
diagrams required to describe exclusive scatterings at tree level. It has the additional 
advantage that the real and imaginary part of eq. ([!]) are generated with the same 
algorithm and that, in principle, the integrand once generated is not manipulated any 
further, hence no mistakes are stirred in. The method consists in evaluating numerically 
the real and imaginary part of the integrals 

f{x) dx (10) 



o x — a + ie 



for a number of small but finite values of e, since for too small values of e numerical 
instabilities take over, and then to extrapolate to e = 0. In || Farrar and Zhang noted 
that for 5. x < e < 1. x 10~ 2 the results seemed to depend very little on e and 
therefore assumed that the amplitudes obtained for e = 5. x 1CT 3 could be taken as a 
good approximation to the e = results. This procedure was later criticized in since 
it failed to reproduce the correct answer in a number of cases in which an analytical 
evaluation was feasible. 

Our findings, for the much simpler set of diagrams we have been studying, are in 
agreement with those in @. For e = 5. x 10~ 3 the difference between the amplitudes 
evaluated with the ie method and with the two others can be as large as 4 5% 
and 20 -r- 30% for the real and imaginary part of the amplitude respectively. For 
e = 5. x 10~ 4 , a value for which the integrals show no sign of numerical problems, 
the difference can still be of several percent. In most much better agreement, 

compatible with the numerical errors, can be obtained fitting a number of results in 
the range 5. x 10 -4 < e < 5. x 10~ 3 with a quadratic polynomial in e whose coefficients 
are determined by the least squares method, and assuming the value of the polynomial 
for e = as the true result. Some caution is however needed. We have found instances, 
for perfectly reasonable values of the scattering angle and of the photon virtuality, in 
which the result is disturbingly sensitive to the chosen range. In some cases, in order 
to reproduce the correct result, the fit has to be performed on a different range, for 
instance 1. x 10 -4 < e < 1. x 10~ 3 . The difference with the result obtained through 
a fit in the range 5. x 10~ 4 < e < 5. x 10~ 3 can be as large as 10%. This obviously 
casts some doubts on the possibility of a completely automatic implementation of the 
method. 

As already noted [jTTJ, the limit in which the mass q 2 = —Q 2 < of the virtual 



photon vanishes is particularly tricky. For non-vanishing photon mass a larger number 
of propagators can go on mass-shell and, after partial fractioning, one discovers that 
different terms require different ranges in e to give an acceptable extrapolation to e = 0. 



Results 

We define the scattering kinematics for e(k, h e )n(p) — ► e(k' , h e )7r(p') r y(q' , A) as in [||. 
The initial and final pion and the final photon have momenta: 

p = P(l, 0,0,-1) 

p> = Q'(l,-sin0,O,-cos0) (11) 
q' = 0^(1, sin 0,0, cos 0) 

while momenta of the initial and final electrons are 

= K(l, cos sin a, sin sin a, cos a) (12) 
fc' = K'(l, cos sin a', sin 0sina', cos a') 

The diagrams contributing to e7r scattering can be divided into two sets. The first set 
includes all diagrams in which the final photon is emitted from the quark lines and 
will be referred to as the pion Compton set. In this case for the momentum of the 
intermediate photon we have 



g=(^P 2 -Q 2 ,0,0,P) (13) 

The second set includes all diagrams in which the final photon is emitted from the 
electron line. 

Momentum conservation requires the following relationships among the parameters 

K sin a = K' sin a! 

P = K cos a — K' cos a (14) 
P + K = 2Q' + K' 
Q 2 = AQ'(P-Q') 

For a (virtual) photon propagating along the positive z-axis, the polarization vectors 
are defined in the standard way 



eZ, R = 4|(0, -1, ±», 0) el = A={$±*$) = 0, 0, ±1) (15) 



Writing the amplitude A^ e for the pion Compton diagrams as: 

< = M v gT W? 2 (16) 

where M v describes the emission of a virtual photon from the lepton line and = 
V tt0 e p * describes the scattering of this photon off the pion, and using the identity 

9 — e R e R + e L e L e + e - ~ e - e + \ LI ) 

one can write 

- q 2 A = M L T L + M R T R -M^T + - M + T_ (18) 



Gauge invariance of the amplitude T M requires 

q ■ T = Jp 2 - Q 2 T° - PT 3 = Jp 2 - Q 2 e ■ T + Pe 3 ■ T 



(19) 



^ /p 2 - q 2 (t+ + r_) + -^p (t+ - r_) 



Vn U ef€"* for incoming photons 



It is therefore sufficient to give the amplitudes 
with polarization e^, e R and e+. 

All the computations have been made in the formalism of [12|. Our results for 
the unintegrated amplitudes coincide with those in ||. In all formulae we use the 
abbreviations c = cos6*/2, c a = cos a/2, c a i = cosa'/2, s = sin 0/2, s a = sin a/2, 
s a > = sina'/2. 

It is amusing to note that the amplitude for rc-f^ - ► t^Ir with real photons has a 
very simple form 

V LR = C (e l -e 2 ) 2 (x-y)s 2 /c 2 (20) 

where Co is a constant, a; and y are the momentum fractions of the initial and final 
quark and e\ and e% are the charges of the two quark lines. This is in agreement with 
a general theorem on helicity amplitudes [O . 



We present only the amplitudes for right-handed final photons, since all others are 
related by parity The helicity amplitudes A^ e for the pion Compton contribution to 
en — ► e7T7, with the definitions given in (|15|), are 



A R 

A* 



c a >c a V +R + s a 'S a V—ji - e^SatCaVnR - e l<t> c a <s a V LR 
c a >c a V +R + s a /s a V- R - e l<p c a is a V RR - e~ t4> s a >c Q V LR 



times the common factor 



ca 



2e^[2KK' 
~~Q 2 



(21) 
(22) 

(23) 



The amplitudes B\ e for the diagrams in which the photon is emitted from the electron 
line are 



B 1 



B 



R 



K'(e i(p s a ,s + c a ,c)(e- i(t, s a ,c - c a is)/2k' ■ q' 

x {Pc a >c a + Q'{e l ^s a is + c a ic){e" l ' l> s a s + c a c 
+ (e~ i(t, s a c - c a s)/(-2k ■ q) x [Pc a ,(Kc a (e^s a s + c a c) - Q'c) (24) 

+Q'(e^s a/ s + c a ,c) (K(e i4, s a s + c a c)(e~^s a s + c a c) - Q' 
K(e i(t) s a s + c a c)(e~^s a c - c a s)/(-2k ■ q') 

X (^Pc a 'C a + Q'fe^SafS + c a rc)(e l ^s a s + c a c)) 
+ (e- i(t, s a ,c - c a ,s)/2k' ■ q' x [Pc a (K' 'c a ,{e i4 's a ,s + c^c) + Q'c) (25) 



+Q'(e icl> s a s + c a c){K' (e^s^s + c a /c)(e~^s Q /s + c Q /c) + Q' 
times the common factor 

2 4 e 3 V2KK J 
Cb — ; P?r ytp P ' ) 



t 



(26) 



where t pp i = —2p-p f and F n is the pion form factor. The full amplitude is simply given 
by + . Combining our results and the data for the pion electromagnetic form 
factor a prediction for the high energy reaction en —>■ en'j can be obtained. 



We have tried three different forms for the pion wave-function denoted CZ [14], P2 



and P3 @: 



<Pcz{x) = 30/^(1 - x ){2x-lf 

<p P2 (x) = 6f„x(l -x) [-0.1821 + 5.91(2x- l) 2 ] (27) 
P3 (x) = 6f„x(l - x) [0.6016 - 4.659(2x - l) 2 + 15.52(2x - l) z 

where f n = .133 GeV. For the strong coupling constant we have adopted the traditional 
value a s = .3. The expressions (|27j) for the wave function have all been derived with 



the sum rule method 16, 17 



The amplitudes for — > -k^r (RR), t^Il — ► t^Jr (LR) and — > (+i?) are 
given in fig.l through 3 for different values of the photon virtuality 77 = Q 2 /4Q' 2 . All 
other amplitudes can be obtained through parity transformations. The corresponding 
phases are given in fig. 4 through 6. 

The amplitudes show a complex dependence on the scattering angle and from the 
photon virtuality. The +R amplitude grows as expected with r\ and only for r\ > .5 
gives a significant contribution. It peaks in the forward and backward direction and 
otherwise has a smooth behaviour with a phase close to 180°. The LR amplitude is 
large in the backward direction and small for small scattering angles as a consequence 
of angular momentum conservation. It has a complicated dependence on 7], decreasing 
for a while for increasing values of rj and then reversing the trend. The RR amplitude 
is the only one with a non trivial phase for a real incoming photon and gives the largest 
contribution. It monotonically increases from small to large scattering angle for real 
photons. For virtual photons it presents a minimum at cos# ~ —0.3, which is quite 
deep for the CZ wave-function, and becomes much larger than the amplitude for real 
photons in the forward direction. Fig. 4 shows that the minimum corresponds to a zero 
of the real part of the amplitude. The RR contribution decreases rather smoothly with 
increasing rj. 

Our results for the LR amplitude and for the real part of the RR amplitude agree 
with the results of [[|, which neglects all imaginary contributions. 

The predictions do not depend too strongly on the choice of the wavefunction within 
the limited set which has been considered. With the exception of the RR amplitude for 
the CZ wave-function and large virtualities, the amplitudes squared change at most by 
a factor of two or three for different wave-functions, in agreement with past experiences. 
The phases are also only moderately dependent upon the wave-function. 



Conclusions 

We have computed pion Compton scattering in perturbative QCD for real and space- 
like initial photons. 

We have studied three different methods for the convolution of the hard amplitude 
with the pion wave-functions. All procedures give the same results when properly 



applied. We find that a careful extrapolation to e = is needed in the ie method. This 
may partially explain some disagreement between calculations performed in the past. 
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Figure Captions 

fig.l. s 3 d ^ s0 for — > 7T7^j in the 717 center-of-mass for the three different wave- 
functions (a:CZ, b:P2, c:P3) and for different values of the photon virtuality 

fig. 2. Same as in fig.l for it^l — ► t^Ir- 
fig. 3. Same as in fig.l for Tvy + — > 777^. 

fig. 4. Phase in degrees for 777^ — > 717^ as a function of the scattering angle in the 717 
center-of-mass for the three different wave-functions (a:CZ, b:P2, c:P3) and for 
different values of the photon virtuality. 

fig. 5. Same as in fig. 4 for 717^ — > 717^. 

fig. 6. Same as in fig. 4 for 777+ — > 777^. 



